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ABSTRACT 

AGN feedback is believed to play an important role in shaping a variety of observed galaxy 
properties, as well as the evolution of their stellar masses and star formation rates. In particu- 
lar, in the current theoretical paradigm of galaxy formation, AGN feedback is believed to play 
a crucial role in regulating star formation activity in galaxies residing in relatively massive 
haloes, at low redshift. Only in recent years, however, has detailed statistical information on 
the dependence of galaxy activity on stellar mass M*, parent halo Mtjm mass and hierarchy 
(i.e. centrals or satellites) become available. In this paper, we compare the fractions of galaxies 
belonging to different activity classes (star-forming, AGN and radio active) with predictions 
from four different and independently developed semi-analytical models. We adopt empirical 
relations to convert physical properties into observables {Ha emission lines, OIII line strength 
and radio power). We demonstrate that all models used in this study reproduce the observed 
distributions of galaxies as a function of stellar mass and halo mass: star forming galaxies and 
the strongest radio sources are preferentially associated with low-mass and high-mass galax- 
ies/haloes, respectively. However, model predictions differ from observational measurements 
in many respects. All models used in our study predict that almost every Mdm > 10 12 M Q 
dark matter halo and/or A/* > lO 11 ^/© galaxy should host a bright radio source, while only 
a small (few per cent) fraction of galaxies belong to this class in the data. In addition, radio 
brightness is expected to depend strongly on the mass of the parent halo mass in the models, 
while strong and weak radio galaxies are found in similar environments in data. Our results 
highlight that the distribution of AGN activity as a function of stellar mass provides one of 
the most promising discriminants between different gas accretion schemes. 
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1 INTRODUCTION 

In the last decade, a consensus has grown that there is 
a strong correlation between the presence of active galac- 
tic nuclei (AGN) powered by gas accretion onto Super Mas- 
sive Black Holes (SMBHs), and the properties of their host 
galaxies. From the observational viewpoint, well defined cor- 
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spheroidal component of the host galaxy point towards a joint evo- 
lution of these two components and a self regulated growth of 
the system. The details of the relative growth of the two compo- 
nen ts are, however, still uncertain (fo r a critical discussion, see 
e.g. lPeng et alj2006l ; IZheng et al.ll2009t) . 

From the theoretical viewpoint, AGNs have long been thought 
to represent a transition phase in galaxy evolution, with no sig- 
nificant influence on the evolution of the physical properties of 
galaxies as a function of cosmic time. Nowadays, this view has 
changed, and AGNs are viewed as a crucial "ingredient" in galaxy 
evolution. Very luminous AGNs (i.e. quasars) release tremendous 
amounts of energy on very short time-scales, eventually trigger- 
ing powerful galactic winds, that can deplete the cold gas reser- 
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voir of the host galaxies, and quench their star formation activity 
jPi Matteo et aljl2005l: iMonaco & Fontanotll200^ ; IPounds & Page! 
120061 : iGanguIv & Brothertonll2008f) . In addition, in order to re- 
produce the red colour and low levels of star formation of local 
massive galaxies, a viable physical mechanism is required so as 
to quench the predicted strong cooling flows in relatively massive 
haloes at low redshift. Low lum inosity radio AGNs provide an 
elegant solution to this p roblem dBower et alj|200fj ; ICroton et al .1 
120061 : Isii acki et al"]|2007l) : radio jets, accelerated by the central en- 
gine, can efficiently transport energy from c entral regions outward , 
offsetting cooling onto the central regions 1 Binnev & Taboj[T995l ; 
iRuszkowski & Begelmanll2002l ; lKaiser & Binney||2003l) . 

In galaxy formation models, high accretion events are usu- 
ally associated with galaxy mergers, that can destabilise large 
amounts of cold gas in the colliding galaxies and drive it to- 
wards the centre of the merger remnant. Observational data indi- 
cate a strong correlation between mergers and e nhanced star forma- 
tion rates in ultra luminous infrared galaxies 1 Sanders & Mirabell 
19961: iBarton et alj|200(j|; IPasquali et"ai]|2004lWoods et alj|2006l : 



Lin et alj|2007l : iLi et alj|2008r) . At the same time, theoretical stud 



ies have pointed out that the rapid inflow of cold gas towards 
the centre may feed the central SMBH and produce a luminous 
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quasa r phase ( Negroponte & White 
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1994 iMonaco et alj |2000| : iGranato et al.l T2OOII : ICattaneo et al 
2003) . The feedback arising from this process is thermal, i.e. origi 



nating from the coupling of a fraction of the QSO bolometric lumi- 
nosity with the interstellar medium. In addition, simulations show 
that, in the case of a merger between two almost equal-mass galax- 
ies, the re mnant morphology is spheroidal ( i.e. dominated by ran- 
dom orbits. lBarnesll992llWalker et al.ll996T ), Therefore, if the bulk 
of SMBH growth and the formation of the galactic spheroids are 
both linked to major mergers, this scenario provides an elegant ex- 
planation for the str ong correlations observed between these two 
galac t ic components |Kauffmann & Haehneltl2000l ; lvolonteri et al .1 
l2003l : lHopkins et al.ll2005l,l2006h . 

Another form of feedback, associated with radio activity, is 
believed to play an important role in quenching cooling flows at 
the centre of massive clusters. The development of radio jets might 
be connected with low accretion rates, analogo usly to what hap- 
pens for X-ray binaries (see e.g. Ijestea 120050 . In this scenario, 
the dominant mode for AGN feedback is the injection of me- 
chanical energy into the ICM, associated with the development of 
X-ray cavities inflated by t he relativistic jets jBirzan et al1l2004 
iMcNamara & Nulsenl l2007h . The physical state of the accretion 
flow, as well as the spin of the SMBHs are thoug ht to play a role 
in det ermining the accretion rate efficiency (see e.g. lFanidakis et all 
l2009h . but details about the physical mechanisms leading to the loss 
of angular momentum of the gas, the subsequent accretion of mate- 
rial onto the SMBH, and the conversion of mechanical into thermal 
energy are still to be understood. 

The physical processes mentioned above span from the Mpc 
scale of galaxy mergers to the sub-pc scale of accretion discs and 
their dependence on physical properties of galaxies (i.e. gas frac- 
tions, gas density, morphological type) are still uncertain. There- 
fore, a semi-analytic approach that couples a statistical (or numer- 
ical) description of the growth of the cosmic structure with sim- 
ple (yet physically and/or observationally motivated) prescriptions 
for the various physical processes at play, provides an efficient and 
valid tool to test specific physical assumptions and to efficiently 
explore the parameter space defined by each model. 

The recent ana lysis of the SDSS data carried out by 
IPasquali et a 1. (2009) has opened new interesting prospects for 



a detailed comparison with predictions from theoretical models. 
This ana lysis has taken ad vantage of the group catalogue de- 
fined by I Yang et alj d2007l) applying a halo-based group finder 
dYang et aL 2005 ) to the New York University Value-Added 
Galaxy Catalogue dBlanton et all 120051) . based on SDSS-DR4 
dAdelman-McCarthv et alj 120060 . The group finder provides an 
estimate of the parent Dark Matter (DM) halo mass (Mum) of 
each galaxy group/cluster, as well as of the "hierarchy" (i.e. the 
central or satellite nature) for all galaxies. Stellar masses esti- 
mates (M^) are computed using the observed relatio n between 
the s tellar mass-to-light ratio and the galaxy colour dBell et al] 
|2003|). The same group ca t alogue has been used in a series of pa- 
per dWeinmann et alj2006llvan den Bosch et alj2008 l: lKimm et alj 
l2009l ; IWeinmann et al.l2009l ; |Pasquali et alj20ld) to investigate the 
dependence of galaxy pr operties on s t ellar m ass, environment and 
hierarchy. In particular, iKimm et alj d20090 have used the same 
group catalogue to study the empirical dependence of star forma- 
tion activity on internal galaxy properties and large-scale environ- 
ment, and have compared observational measurements with results 
from semi-analytic models. These studies have shown that most 
models including some prescriptions for AGN feedback reproduce 
reasonably well the observed dependence of the red/passive frac- 
tion as a function of stellar mass and halo mass for central galaxies. 
The same models, however, produce too large a fraction of red and 
passive satellite galaxies. 



IPasquali et aL I d2009l P09 hereafter) have extended this analy- 
sis by matching the group catalogue with th e catalogue of emis- 
sion line fluxes by iKauffmann et alj d2003l) . with the National 
Radio Astronomy Observa tory Very Large Array Sky Survey 
(NVSS: ICoiidon et al.|[l998h . and th e Faint Images of the Radio 
Sky at Twenty centimetres (FIRST; iBecker et al. I [19951) . P09 de- 
fined four different "activity" classes. Star-forming galaxies (SFGs) 
were separated from AGN host galaxies using line-flux ratios 
[Oml A5007/ff/3 and [jVii 1A6584/ffa (the so-called BPT dia- 
gram, [BaldwirretaL] [l98l]). Radio sources detected at 1.4GHz 
were split into sources with low (LRA) and high (HRA) radio ac- 
tivity. Using this matched catalogue, P09 studied the environmen- 
tal dependence of galaxy activity, as well as its dependence on the 
galaxy stellar mass. They showed, in particular, that SFGs are pref- 
erentially low-mass galaxies living in relatively small haloes, op- 
tical AGN activity is linked to intermediate mass galaxies and en- 
vironments, while radio galaxies mainly reside in massive cluster 
galaxies. Given the strong stellar-halo mass correlation for central 
galaxies, it is unclear whether the activity levels for these galax- 
ies are more strongly related to their stellar or halo mass. As for 
satellite galaxies, P09 found that the dependence of their galaxy 
activity on halo mass is about four times weaker than that on stellar 
mass, i.e. activity in satellites is suppressed with respect to central 
galaxies of similar mass. 



In this work, we present a detailed comparison between the 
P09 observational data and predictions of four different semi- 
analytic models. This paper is organised as follows: in sec.|2]we de- 
scribe the different semi-analytic models we consider in this study, 
while in sec. [3] we discuss our strategy for converting model pre- 
dictions into observational quantities. In sec. [4] we then compare 
these predictions with the results presented in P09. Finally, we dis- 
cuss our results and give our conclusions in sec. [5] 
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Model 


n 


n A 


h 


(78 


Pasauali et al. (2009) (P09) 


0.238 


0.762 






MORGANA 


0.24 


0.76 


0.73 


0.80 


Kane et al. (2005) 


0.25 


0.75 


0.71 


0.75 


Somerville et al. (2008) 


0.28 


0.72 


0.70 


0.817 


Wans et al. (2008) 


0.226 


0.774 


0.743 


0.722 



Table 1. Cosmological parameters adopted by the semi-analytic realiza- 
tions. 



2 MODELS 

In this paper, we consider predictions from four independently 
developed theoretical models for galaxy formation. In particu- 
lar, we consider the most recent implementat ions of four differ- 
ent Semi-Analytic Models (SAMs, see e.g. iBaughl 120061 for a 
recent review) for galaxy formation and e volution, in the con - 
cor dance ACDM cosmology: MO RGANA dMonaco et alj 120071) . 
the iKang & van den Bosch] d2008l) model (KVB08 he reafter), the 
ISomerville et alj d2008l) model (SHC08 hereafter), the I Wang et all 
d2008h model (WDL08 hereafter). In these models, the evolution 
of the baryonic component is followed by means of an approx- 
imate description of the physical processes at play (such as gas 
cooling, star formation and feedback) and of their interplay with 
gravitational processes (i.e. dynamical friction, tidal stripping and 
two body mergers), linked to the assembly of the large scale struc- 
ture of the Universe. These "recipes" include a number of parame- 
ters which are usually fixed by comparing model predictions with a 
set of low-redshift observations. Despite their simplified approach, 
SAMs have turned into a flexible and widely used tool to explore a 
broad range of specific physical assumptions, as well as the inter- 
play between different physical processes. 

We refer to the original papers for a detailed description of 
the modelling adopted for the vari ous physical processes consid- 
ered (see also iFontanot et alll2009l for a comparison between the 
MORGANA, WDL08 and SHC08 models). In the following, we 
will focus on the modelling of gas accretion onto SMBHs and 
different modes of AGN feedback. All models are based on an 
analytical or numerical description for the redshift evolution of 
the mass and number density of dark matter haloes (i.e merger 
treefl). The cosmological parameters adopted are slightly differ- 
ent (see table [T); we make no attempt to correct the results, since 
we expect these s mall differences to have negligible influence on 
model predictions dWang et"al]|2008h . In all cases, the mass reso- 
lution is sufficient to resolve galaxies with stellar mass larger than 
10 9 Mq. In the following, we parameterise the Hubble constant as 
H = lOO/tkms^Mpc -1 . 

Gas inside DM haloes cools and condenses primarily via 
atomic cooling and thermal Bremsstrahlung, and forms a rota- 
tionally supported disc. Star formation is modelled using simple 
empirical (Schmidt-Kennicutt-like) recipes; supernovae and stellar 
winds deposit both thermal and kinetic energy into the cold gas and 
might re-heat or expel part of this gas; galaxy mergers might trigger 
"bursts" of star formation and play an important role in triggering 



The WDL 08 and KVB08 mo dels use merger trees ex tracted from N-body 
simulations jWang et al .|l20()8| and ljing & Sutcll2002l respe ctively); MOR- 
GANA uses the Lagrangian semi-analytic code PINOCCHIO dMonaco et alj 
|2002|) . and SHC08 use a method based on the Extended Press-Schechter 
formalism, described in ISomerville & Kolattl fl999l> . 



AGN activity. The details of the modelling of these processes differ 
among the four SAMs we consider in this paper. 

We recall that the WDL08 and KVB08 models explicitly fol- 
low the evolution of DM substructures, until tidal truncation and 
stripping reduce their mass below the resolution l imit of the sim- 
ulatio n (for details, see |Pe Lucia & Blaizotll2007l and IKang et alj 
120051 respectively). Satellite galaxies are then assigned a residual 
survival time, which is roughly twic^l the value computed us- 
ing th e classical dynamical friction formula dBinnev & Tremainel 
Il987l) . SHC08 and MORGANA model the evolution of DM sub- 
structure s using analytic recipes based on the fitting formulae pro- 
posed by |Boyi^Mchinet^|200i) and lTaffoni etai1d2003l) . re- 
spectively. A de tailed comparison betw een these different formula- 
tions is given in lDe Lucia et alj J2010h . 

All models follow the growth of SMBHs at the centre of model 
galaxies, and differentiate between the so-called "bright-mode" (or 
"QSO-mode") and "radio-mode" associated with the efficient pro- 
duction of radio jets. In the following, we discuss in more detail the 
modelling adopted for these physical processes. 



> IWangetalJd2008h . This model represents a generalisation of 
the |Pe Lucia & Blaizotl d2007l) mo del to the WMAP3 cosmology 



IWang 



et al 



( 2008 1) simulations. P redictions from 



assumed in the 

the model presented in iDe Lucia & Blaizotl d2007l) applied to the 
Millennium Simulation are publicly availabltQln this model, the 
"bright-mode" is triggered by gala xy-galaxy mergers following th e 
approach originally proposed by Kauffmann & Haehneltl d2000h . 
and extended to minor mergers bv lCroton et al. 1 20061) : during each 
merger (both minor and major) a fraction /bh of the cold gas 
M co id in the progenitors is assumed to be accreted on the SMBH 
at the centre of the merger remnant (the pre-existing SMBHs are 
assumed to merge instantaneously). The total accreted mass during 
the merger event is parameterised as: 



A„,C06 r /- 

Am QM — /BH C,m 



1 + 



l 280km/s J 



(1) 



where £ m represents the baryonic ratio between the progenitor 
galaxies (defined such as £ m < 1). According to the above equa- 
tion, equal mass mergers result in stronger accretion rates than mi- 
nor mergers between galaxies with similar amounts of cold gas. 
This is to be expected as numerical simulations show that major 
mergers are responsible for a stronger dynamical response and a 
more rapid infall of gas towards the centre, /bh = 0.03 is a 
free parameter, chosen by requiring the model to fit the local re- 
lation between the B H mass and the mass of the hosting spheroid 
dHaring & Rixl2004) . 

"Ra dio-mode" feedback is implemented as in ICroton et alj 
d2006h and is assumed to result from the accretion onto the central 
BH of hot gas from a static atmosphere. The accretion is assumed 
to be continuous and is described by the phenomenological recipe: 



• C06 

m RM 



^C06 



Mb 



/hot \ / V v i 



10 8 M Q y V 0.1 / ^OOkm/s 



(2) 



where /hot is the fraction of the total halo mass in form of hot 
gas and K poe = 6.5 x 10~ 6 Mo/yr. We note that in their origi- 
nal paper, ICroton et al] d2006l) demonstrate that this formulation is 



2 A fudge factor of ~ 2 accounts for recent findings that the classical for- 
mula systematically under-estimates merger times computed from numeri- 
cal simulations (e.g. lBovlan-Kolchin et al, 2008; Jiang et al. 2008) 

3 see http://www.mpa-garching.mpg.de/millennium/ 
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consistent with e xpectations of theore tical models based on the an- 
alytical results of lBertschingeJjT989l ) and on a Bondi-Hoyle model 
for BH feeding. The amount of energy injected by the AGN into the 
intra cluster medium (ICM) through this channel is: 



r C06 _ ■ C06 2 
L BM — ™RM c 



(3) 



where e = 0.1 represents the conversion efficiency of rest mass 
to energy. The efficiency of gas heating by AGN is assumed to be 

inde pendent of the halo mass. 

• iKang & van den Boschl i2008t). The "bright-mode" is mod- 



Kauff mann & Haehnel and is associated with 



elled as in 
major and minor mergers: 

M co i d 



A K06 

Am-QM = ffBH 



1 + 



280km/; 



:)" 



(4) 



g BH = 0.03 is a free parameter, again calibrated using the ob- 
served relation between the mass of the BH and the mass of hosting 
spheroid. 

The "radio-mode" is modelled assuming that the energy injected 
from the central SMBH is proportional to its Eddington Luminosity 
L e dd and that the heating efficiency scales as a power law of V v ir 
(Kang etal. 200j]): 



rK06 



ftK06 L c 



(~20C 



(5) 



where kko6 ~ 2 x 10~ j is a free parameter, fixed by requiring the 
model to fit the bright end of the local galaxy luminosity function. 

Another important feature of this model is the treatment of 
hot gas associated to recently accreted satellites. The other mod- 
els adopted in this study assume that the hot gas reservoir is in- 
stantaneously stripped when galaxies are accreted onto larger sys- 
tems and become sate llites. As argued by many authors (see e.g 
Weinmann et al .l2006h . this assumption is probably too strong, and 
results in red satellite fractions significantly higher than observed. 
KVB08 show that by assuming a constant stripping rate over a 
timescale of ~ 3 Gyrs, the model is able to predict a fraction of 
blue sate llites which is in qualitative agreement with observations 
(see alsolWang et al J2007h. 

• ISomerville et alj J200S|) . The model used in this study is dis- 
cussed in detail in SHC08, but u s es the updated bulge formation 
model discussed in lHopkins et al. I J2009h . Other minor changes in- 
clude a new calibration of the galaxy formation parameters, ap- 
propriate for a WMAP5 cosmology. In this model, the "bright- 
mode" is associa ted with merger even ts, following the prescrip- 
tions proposed bv lHopkins et al. I J2007I) and based on a suite of N- 
body simulation s of galaxy-galaxy mergers dRobertson et al .120061 : 
ICox et al I l2006h . For each merger event, the ma ss of the final 
sphero idal component M sp h is computed following iHopkins et al.1 
( l2009h . and the final mass of the SMBH is estimated using the 
simulation-calibrated relation dHopkins et alj|2007h : 



log(M BH /M B ph) 



-3.27 + 0.36 erf [(/ gas - 0.4)/0.28] (6) 



where / gas represents the "effective" gas fraction (defined as the 
sum of the cold gas masses in the merging galaxies, divided by 
their total baryonic masses). The accretion of gas onto the SMBH 
is modelled throug h the generalised func tional form for QSO light 
curves proposed bv lHopkins et al .1 j2005l) : the SMBH grows at the 
Eddington rate, until it reaches a critical mass M B H,crit. Further 
accretion is then reduced by the onset of a pressure driven outflow 
(which eventually halts the accretion) and this regime is described 
by a power-law decline. The critical mass is defined as: 



M B H,crit = 1.07/ BH ,crit (M BH /lO 9 Af ) 



(7) 



where / B H,crit = 0.4 is a free parameter. In this scenario, bright 
QSOs spend most of their lifetime in the power-law decline phase, 
but most of the BH mass is accreted at the Eddington rate. If the 
mass of the SMBH is already larger than the critical mass, no accre- 
tion is allowed on the central object. This might happen frequently 
at low redshift, where the most massive spheroidal galaxies have 
already assembled most of their stellar mass. 

The "radio-mode" feedback implementation is based on the as- 
sumption that low-accretion events are fuelled by Bondi-Hoyle ac- 
cretion: 



Wl bondi = 7T (GAf BH ) p(r\) C s 



(8) 



where c s represents the gas sound speed and p{rx) is the density 
of the hot gas at the accretion radi us. By combining eq. [8] with 
the isothermal cooling flow solution dNulsen & Fabianll200of) . it is 
possible to rewrite the accretion rate as: 



■ S08 
m RM 



KS08 



kT 



I- 



Mb 



A(T,Z h ) J \l0 s M t 



(9) 



where T is the virial temperature of the hot gas, k is the Boltzmann 
constant, and A(T, Zh) represents the cooling rate, ksos = 3.5 x 
10~ 3 is a free parameter. The energy that is effectively injected into 
the ICM (the heating rate) is obtained by converting the accreted 
matter into energy: 



r S08 _ ^,S08 2 



(10) 



In the standard model, r% eB t = 1.0 and e = 0.1. 

• MORGANA. In this mode l, the grow t h of S MB Hs is modelled 
follow ing t he prescriptions by Umemural d200lh and lGranato et al.1 
J2004I) (see lFontanot et alj|2006f for more details on the actual im- 
plementations of these prescriptions in MORGANA). The model 
follows the evolution of a cold gas reservoir around the SMBH, 
which is assumed to form through the same physical processes that 
drive the formation of the bulge. The formation rate of the reservoir 
m resv is then assumed to be proportional to the star formation rate 
in the bulge rh*,B: 



— fresv 1^1* 



I 100M G 



(11) 



where / rcS v = 3 x 10~ 2 and a Icav — 1 are free parameter^, 
fixed by requiring the model to re produce the < z < 5 evolution 
of the QSO luminosity function dFontanot et al.l [2006h . The cold 
gas mass in the reservoir Af rcsv can be accreted onto the SMBHs, 
once it loses its residual angular momentum, through turbulence, 
magnetic field or radiation drag. The BH accretion rate is assumed 
to be determined by the viscosity of the accretion disc and is written 
as dGranato et aL1l2004) : 



• M07 _ 

m BH — ^visc 



G 



M B 



3/2 



1 + 



Mb 



1/2 



where as is the velocity dispersion of the bulge, and k v 



(12) 



10" 



is cho sen following the theoretical arguments by iGranato et al.l 
d200lh . The accretion rate is limited at the Eddington rate. In this 
model, it is the rate of accretion that determines the 'nature' of 
the feedback: the "bright-mode" is associated with accretion rates 



4 Note that with this parameter choice, the equation simplifies to: m re av 

/resv B- 
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larger than 1 per cent of the Eddington limit, while the "radio- 
mode" is sustained by lower accretion rates. For the radio-mode, 
MORGANA assumes a bolometric luminosity: 

r M07 ^,M07 2 

The efficiency (kmo7) at which the emitted energy heats the hot 
halo component is assumed to scale with the circular velocity of 
the host halo and, in the standard implementation, is normalised to 
the value at 10 3 km/s 

KMm = ( lOOOWs ) (14) 

For consistency with the other models, in the following we neglect 
this virial velocity dependence (i.e. we set kmo7 = 1-0). We have 
verified, however, that its inclusion does not change significantly 
the results discussed below. 

For high accretion rates, the model assumes that only 10 per cent 
of the bolometric luminosity is available to heat the halo gas (con- 
sidering that about 10 per cent of quasars are radio loud). We note 
that, at variance with the other models used in this study, in MOR- 
GANA the gas accreted during the "radio-mode" is cold. In addi- 
tion, gas accretion is always associated with some star formation, 
responsible for the final loss of angular momentum of the accreted 
material. As a result, massive central galaxies in this model are too 
blue (because they form too many stars) with respect to the obser- 
vational determinations. 

In this paper, we carry out a comparison between predictions 
from the 4 models introduced above and the observational measure- 
ments presented in P09. We stress that none of the models used has 
been tuned to reproduce the observational data considered. There- 
fore, model results discussed below should be considered as gen- 
uine model predictions. 

This section illustrates that the implementations of SMBH 
growth and AGN feedback adopted in the models considered in 
this study, differ in a number of details. We note, in particular, that 
minor mergers do trigger accretion onto BHs in three models used 
in this study (MORGANA, SHC08 and WDL08), satellite-satellite 
merger are considered in KVB08 and WDL08, while only MOR- 
GANA considers also accretion of gas triggered by disk instabili- 
tiefl In addition, AGN driven winds are explicitly included only in 
MORGANA and SHC08. The SHC08 model includes a more sophis- 
ticated treatment for the evolution of the AGN during the "bright" 
phase, explicitly including a description for the AGN light curve. 

The models also share a number of similarities. "Bright- 
mode" feedback is tightly linked to galaxy-galaxy mergers in 
WDL08, KVB08 and SHC08. MORGANA adopts a different 
scheme, but also for this model the high accretion rates needed in 
the "bright-mode" regime are often associated with galaxy-galaxy 
mergers. In addition, three of the four models used in this study 
(WDL08, KVB08 and SHC08) assume (either directly or indirectly 
as in the SHC08 model) a strong relation between the mass of the 
parent DM halo and the strength of the "radio-mode" feedback. 

The variety of assumptions made to model these processes re- 
flect the fact that they are still not well understood, and that almost 
all details of the physical mechanisms acting on sub-pc scale (i.e. 
the final stages of angular momentum loss, the properties of the in- 



5 We note th at disk instabilities represent the main channel for SMBH 
growth in the Bower et a]] J2006h model, which is not considered in this 
study. 
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Figure 1. Selection criteria for the AGN activity class. Black solid contours 
refers to the 5%, 10%, 30%, 50%, 70%, 90% number density levels corre- 
sponding to the observed AGNs in the P09 sample. Yellow, magenta, green, 
blue dots refer to the predictions of the MORGANA, KVB08, WDL08 and 
SHC08 models respectively. Blue dashed line represents the observed flux 
limit converted to intrinsic luminosity at z = 0.1; blue solid line corre- 
spond to the S/N cut of the data. The shaded area corresponds to the final 
selected region. 

ner accretion disc, the AGN light curves, and the collimation of the 
radio jet) are still largely unconstrained. 



3 CONVERTING MODEL PREDICTIONS INTO 
OBSERVED QUANTITIES 

In order to compare model predictions with the P09 results, we 
need to convert the predicted accretion and star formation rates 
into 'observables'. These include, in particular, the luminosity of 
the Ha and OIII lines, and the radio power at 1.4 GHz Pi.4GHz 
used for the estimate of the star formation, AGN and radio activity, 
respectively. All model predictions presented below are constructed 
by considering a cosmological box centred atz — 0.1, and convert- 
ing predicted intrinsic luminosities into apparent fluxes. 

3.1 OIII line luminosity 

Predicted accretion rates onto SMBHs are converted into an esti- 
mate for the luminosity of the OIII 5007A line (Loin) as illustrated 
below. For each model galaxy, we compute the expected bolomet- 
ric luminosity associated with each accretion event in the "bright- 
mode" as: 

ibolo = eiiiQMC 2 (15) 

where we assume e = 0.1. We then estimate Lp m using an em- 
pirical approach based on the results obtained by iHeckman et al.l 
(2005) from a sample of type I and II AGNs. These authors esti- 
mate the mean luminosity ratio between the hard X-ray luminosity 
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(Lux) and Loin for a sample of X-ray selected local AGNs, and 
find that: 



log(L H x/ioin) ~ 2.15 



(16) 



To make use of this result, we first convert bolometric luminosi- 
ties into Lhx by me ans of the bolometric correction proposed by 
iMarconi et al, I J2004l) . We also test an alternative conversion based 
on the bolometric correction proposed bv lWul l l2009l) for a sample 
of type I AGNs: 



log(Lboio) = 0.951og(L O in) + 5.39 



(17) 



and we find that the two approach es provide similar results. We 
note that the iHeckman et alj i 20051) sa mple includes both Type I 
and type II AGNs. while fhe lWul d200^) sample includes only Type 
I AGNs. Since type I AGNs were explicitly excl uded in the P09 
analys is, we will show only results based on the IHeckman et al" 



2005) conversion in the following. In addition, following Simpson 



2005 ), we estimate the Type I fraction as a function of Loin and 



we use it as a statistical correction to account for the removal of 
Type Is in P09. We have verified that this correction does not affect 
our conclusions. 

In order to compare model predictions with the P09 results, 
we also need to apply similar 'selection criteria'. The main require- 
ment for a galaxy to belong to the AGN class in the P09 analysis is 
the detection of the OIII emission line in the SDSS spectra, above 
a well defined signal-to-noise level S/N > 3. This implies that, if 
the OIII flux is relatively weak compared to the continuum, and the 
galaxy is at a relatively large distance (so that its total flux is low), 
the continuum will be noisy and the emission line will not be de- 
tected at sufficient S/N. In fig.[U we plot the distribution of model 
galaxies (coloured symbols) in the Loin-stellar mass plane, and 
compare it with the distribution measured from the SDSS: black 
contours show the 5, 10, 30, 50, 70, and 90 per cent number den- 
sity levels from the observed data. This figure illustrates the conse- 
quences of assuming a S/N limit: OIII lines of comparable lumi- 
nosity are more easily detected in less massive galaxies, due to the 
lower continuum level. Using the distributions shown in fig.Q] we 
define a composite selection criterion for our model galaxies. We 
first convert the observed flu x limit (/qui = 5x 10~ 17 erg/s/cm 2 , 
de-reddened as described in iKauffmann et alj|2003l) into an intrin- 
sic luminosity limit at z — 0.1. This limit is marked with a blue 
dashed line in fig. Q] We then define a cut in specific luminos- 
ity, Loin/M* > 10 5 ' 6 L Q /M Q , which mimics the S/N selec- 
tion (blue solid line) and is chosen to approximately follow the 
5 per cent number density level of the observed distribution. The 
shaded area in fig. Q] highlights the region occupied by selected 
model galaxies. 

We note that predictions from the four models used in this 
study populate different regions of the Loin-stellar mass diagram. 
This suggests that, in principle, we could use these distributions to 
discriminate between competing models. Unfortunately, however, 
observational data do not cover a region of this space that is wide 
enough to provide conclusive evidence in favour of or against a 
given model. As an example, WDL08 and KVB08 show similar 
results in the high luminosity region: this is expected, since these 
two models use very similar prescriptions for gas accretion during 
merger events. In MORGANA, the distribution of AGNs shows a 
strong dependence on M*, similar to that obtained for radio galax- 
ies (see sec. 13. 3\ . We remind that in this model gas accretion is 
treated in the same way in both regimes, the only difference being 
the amount of accreted gas. The correlation between Loin and M* 
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Figure 2. Selection criteria for the SFG activity class. Symbols, lines and 
shading as in fig.[TJ 



originates from the strong dependence of gas accretion on Mbh 
and gb- 



3.2 Star formation rate and Ha line strength 

Predicted star formation rates are converted into Ha line luminosi- 
tie s (Z/Ha) u s ing th e observationally calibrated relation proposed 



bv lKermicuttUl998l) : 

SFR(M e /yi) = 7.9 x 10~ 42 L H a (erg/s) 



(18) 



Since the Ha luminosity provides, on average, an estimate of only 
half the total SF R, assuming an a verage 50% escape fraction of 
ionising photons dKennicutnl 19981) . we reduce the Lh q estimated 
using the above equation by a factor of two. In fig. [2] we show 
the distributions of predicted Lua luminosities, and compare them 
with the observed distribution. In order to compare model predic- 
tions with the data presented in P09, we adopt a strategy similar 
to that illustrated in sec. 13.11 We first convert the apparent flux 
limit (/Ha = 8 x 10~ 16 erq/s/cm 2 de-reddened as described in 
Kauffmann et al. I l2003h . into an intrinsic luminosity at z = 0.1 
(blue dashed line). We then define a cut in specific luminosity 
Lu a /M* > 1O _3 - 9 L /M (blue solid line), that is chosen to 
approximately follow the 5 per cent density level and mimics the 
S/N > 3 cut adopted in the observational analysis. As in fig. \T\ 
the shaded region highlights the selected area. We do not attempt 
to take into account oth er possible biases due to fibre size (see e.g. 
iBrinchmann et al.l2004l for a discussion): we remind the reader that 
the observational estimates are based on fibre spectra, while SAMs 
provide total integrated SFRs for model galaxies. Predictions from 
the four models used in this study populate very similar regions of 
the Lhq - M* diagram. 
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Figure 3. Selection criteria for the Radio activity classes. Symbols, lines 
and shading as in fig.[T] Red dot-dashed line represents the assumed thresh- 
old between LRA and HRA classes. 



3.3 Radio 

In order to estimate the radio emission associated with the "radio- 
mode", we assume that an energy corresponding to Lrm is 
converted into mechanical work that radio AGNs exert on the 
surrounding mediu m. We then use the relation calibrated by 
IWillott et all Jl999h to determine the luminosity at l51Mhz (X151, 
in units of lO^WHz^sr- 1 ): 



L kin = 3 x 10 45 / 3/2 ^!erg/s 



(19) 



In eq.[T9] / is a correction factor that takes into account the system- 
atics in the observat ional measurements, a nd we as sume the aver- 
age v alue / — 15 dHardcastle et al. I [20071 see also lShankar et alj 
120081) . We then compute the 1.4 Ghz luminosity (Li,4Ghz) assum- 
ing a fixed spectral shape /„ oc v ' 7 , and apply the same detection 
limits of P09 (.L 1AQhz > 10 21 h _a WHz _1 ). The full conversion 
between Lrm and L1.4 can be thus expressed by the following ex- 
pression: 



L1.4 = 4tt x 10 28 



7/6 

^3 x 10 45 / 3/2 J X 



10 7 



1.4 
0.151 



0.7 



WHz' 1 



(20) 



We have also tested an alternative me thod based on the empirical 
conversion calibrated on the results bv lBest et alJfeOOd) : 



10 3e W 



= 3 x 



Ll.4GHz 

lO^WHz- 1 



(21) 



and verified that this does not affect our conclusions. 

As in P09, we split radio galaxies into two samples 
with low and high radio activity (LRA and HRA respec- 
tively) using a specific luminosity threshold of Li.4Ghz/M* = 
10 12 h" 2 WHz _1 Lo/M Q . We stress that in the following we do 
not consider radio emission associated with high accretion rates in 



the models, i.e. we n eglect the radio loud fraction of QSOs (see 
e.g. ljiang etZll2007l) . Our estimates of the fraction of radio galax- 
ies are therefore lower limits: however we want to test the con- 
servative hypothesis that all radio activity is linked to the "radio- 
mode" feedback. In fig. [3] we compare model predictions with the 
observed radio power distribution. In order to mimic the obser- 
vational selection criteria, we apply the same apparent radio flux 
limit of 2.1mJy. In fig. [3] we show the corresponding value for 
the intrinsic power at z = 0.1, as a blue dashed line, while the red 
dot-dashed line marks the threshold between the LRA and HRA 
classes. The selection criteria illustrated above influence our re- 
sults, particularly for the LRA sources, but do not affect our con- 
clusions. We will come back to this issue later. We note that models 
predict the correct dependence of PiAGim on M*, with more mas- 
sive galaxies hosting brighter sources, but the correlation between 
radio power and stellar mass appears to be stronger in the models 
with respect to the observational results by P09. The strong de- 
pendence of Pi.4Ghz on Mi, in the WDL08 and KVB08 models 
results from the strong correlation assumed between the strength 
of "radio-mode" feedback and the mass of the parent DM halo. 
Since this mode is active only for central galaxies, and the stel- 
lar mass of these is strongly related to the parent halo dark matter 
mass, this results in a strong correlation between the intensity of 
radio-power and the galaxy stellar mass. Although SHC08 adopt a 
different model governing the efficiency of radio mode activity, as 
already shown in SHC08 their model produces a nearly identical 
dependence on halo mass as that used in WDL08. Therefore we 
expect a similar dependence on stellar mass as well. Among the 
models used in this study, MORGANA shows the weakest depen- 
dence of Pi.4Ghz on M*. This is due to the fact that in this model, 
gas accretion onto SMBHs in the "radio-mode" regime is not di- 
rectly linked to Mom, but is mainly regulated through Mbh and 
m*,s (eq.[72l. 



4 RESULTS 

Following P09, we consider 4 different activity classes: SFG, AGN, 
LRA and HRA (see previous sections), and consider separately 
central and satellite galaxies. For each class, we define the fraction 
f(A\P) of model galaxies belonging to the activity class A and 
with a given property P (i.e. stellar mass and parent halo mass). 
We note that in the P09 analysis, SFGs and AGNs are separated 
using the BPT diagram. In this paper, we simply assume that all 
model galaxies with a detectable OIII line belong to the AGN class. 
Therefore, our model predictions for f(AGN\P) represent an up- 
per limit, since AGNs hosted in SFGs may be misclassified in the 
observations. However, as we will see in the following, we do not 
expect this contamination to affect our conclusions. We note that, 
although the SFG and AGN classes are mutually exclusive (as in 
the real data), the same objects might be classified as radio sources, 
both for model predictions and for observations. 

In fig.|4]and |5] we show the /(A|M D m) and f(A\M*) frac- 
tions, for central and satellite galaxies respectively. Model predic- 
tions are shown using coloured histograms (yellow, red, green, blue 
for the MORGANA, KVB08, WDL08, SHC08 models respectively), 
while shaded histograms show the observational measurements by 
P09. When constructing the model histograms, we apply additional 
cuts to our model samples by considering only galaxies with pre- 
dicted stellar masses between 10 9 ' 2 < M*/i -2 /Af Q < 10 12 and 
living in DM haloes of 10 11 ' 4 < A/dm/i~7 m © < 10 15 , as in 
P09. Typical uncertainties in observed fractions are of the order 
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Figure 4. Fraction of central galaxies belonging to the different activity classes as a function of stellar mass (right panel) and as a function of the mass of the 
parent DM halo (left panel). In all panels the yellow, red, green and blue histograms refer respectively to the MORGANA, KVB08, WDL08 and SHC08 models 
(histograms are slightly offset for clarity). Shaded areas refer to the fractions observed by P09. 
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Figure 5. Fraction of satellite galaxies belonging to the different activity classes as a function of stellar mass (right panel) and as a function of the mass of the 
parent DM halo (left panel). Lines and shading as in fig. [4] 



of 5% for the SFG and AGN classes, 10% for HRA/LRA centrals, 
20% in f{A\M OM ) for HRA/LRA satellites, and 30% in f{A\M*) 
for HRA/LRA satellites. 

We first consider central galaxies (fig.|4]l: the peak of the mock 
distributions shift from low-mass haloes for SFGs to high-mass 
haloes for HRA sources, in agreement with the observational mea- 



surements. In addition, we find that the distributions as a function 
of galaxy stellar mass are quite similar to those computed as a func- 
tion of halo mass. This is expected because of the strong correlation 
between halo mass and s tellar mass for central galaxies in SAMs 
(see e.g. lYang et alj2008h . 

The figure, however, illustrates several discrepancies between 
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model predictions and observations. The KVB08, SHC08 and 
WDL08 models reproduce the observed decrease of the fraction of 
star forming galaxies with increasing stellar and halo mass. MOR- 
GANA predicts higher f(SFG\Mr>yi) than observed, especially at 
intermediate masses. This is due to the inefficient quenching of star 
formation in massive gal axies, leading to a large fraction of active 
massive central galaxies dKimm et all2009l) . Both f(AGN\M OM ) 
and f(AGN\Mi,) are more skewed towards higher masses with 
respect to the observational measurements. The SHC08 model rep- 
resents an exception, with a very low predicted fraction of AGNs, 
and a distribution that peaks at the low-mass end, both as a func- 
tion of stellar and parent halo mass. This is due to the modelling 
adopted for QSO activity during mergers: this model considers a 
critical BH mass, above which no further accretion is allowed onto 
SMBHs. At low redshift, most massive galaxies already host mas- 
sive BHs, so that most mergers at low redshift happen between 
galaxies whose combined BH masses are larger than the critical 
mass adopted in this model and further accretion onto the rem- 
nant BHs is suppressed. In contrast, the other three models used 
in this study still allow gas accretion onto the SMBH of the rem- 
nant galaxy, as long as the progenitors still have some gas left (i.e. 
if the merger is not "dry"). 

The predicted fractions of radio sources as a function of both 
stellar and halo masses are larger than the observed ones by at least 
one order of magnitude: all models predict that more than 80 per 
cent of central galaxies in relatively massive haloes should host a 
strong radio galaxy. The disagreement with observations is partic- 
ularly large for f(H RA\Mdm)'- almost all models predict that all 
Mom > 10 13 Mq/H haloes are associated with a detectable radio 
galaxy. Also the shapes of the /(L7L4| M D m) and f(HRA\MoM) 
distributions predicted from the models used in this study dif- 
fer from those observed (this holds also for f(LRA\M+) and 
/ (H RA\M±))\ bright and faint sources are associated with larger 
and smaller haloes respectively, while the observed distributions for 
these two classes are very similar. The strong correlation between 
radio activity and halo mass is expected in the KVB08, WDL08 and 
SHC08 models, given the assumed proportionality between accre- 
tion rates in the "radio-mode" and the mass of the parent dark mat- 
ter halo. It is interesting that the typical halo mass that separates the 
HRA and LRA classes is similar in these models, despite the dif- 
ferent dependence on circular velocity adopted. We anticipate that 
the precise value of this "transition" mass depends on the details of 
the conversion of heating rates into radio luminosity (we will come 
back to this in sec. [5}. MORGANA is the only model that predicts a 
significant population of low-luminosity radio galaxies for central 
galaxies in massive haloes, but both the predicted fractions and the 
distributions differ from those observed. We stress that the shape 
of the LRA distribution depends on the 'selection criteria' adopted 
for model galaxies (see sec. 13. 3t . A decrease in the assumed flux 
limit translates into larger samples of faint radio sources: this af- 
fects only the low stellar and halo mass tail of the distributions, 
leading to larger f(LRA\P) values. As a consequence, the peak 
of the model distributions moves towards lower masses, increasing 
the discrepancy with the observational measurements. 

We now consider the f(A\MoM) and f(A\M+) fractions for 
satellite galaxies. Most models predict f(SFG\M+) distributions 
that differ from those observed. The distributions are significantly 
affected by the assumed selection criteria at their low-mass end: 
when considering all satellite galaxies above the minimum stel- 
lar mass observed, MORGANA, SHC08 and WDL08 models pre- 
dict an approximately constant fraction of SFGs as a function of 
M*. MORGANA and SHC08 predict large /(SFGI M*) for high- 



mass satellites, while the lowest fractions of SF satellites are given 
by the WDL08 model. Only the KVB08 model predicts the cor- 
rect dependence of the fraction of star forming satellites on stellar 
mass, although the predicted fractions are larger than observed. It 
is interesting that the same model predicts an almost flat distribu- 
tion for /(SFG|A/dm), while the other models are able to catch, 
at least qualitatively, the slight decrease of SFGs with increasing 
halo mass. These results highlight that the modelling adopted in 
KVB08, despite being able to provide the correct fraction of satel- 
lite as a function of galaxy colours, does not fully reproduce the 
observed SF activity of satellites as a function of halo mass. 

The distribution of f(AGN\M*) is too skewed towards high 
M*. The only exception is again the SHC08 model, which shows 
a small peak of activity at M* ~ 10 Mq, but does not repro- 
duce the observed total fractions of satellites with AGN activity. 
f(AGN\Mr>M) is approximately flat in the models, in qualitative 
agreement with data. We stress that in the MORGANA, WDL08 and 
SHC08 models, by construction, no AGN activity is expected in 
satellites, since none of them allows accretion of new cold gas onto 
satellites and satellite-satellite mergers at the same time. Therefore, 
the galaxy-galaxy mergers responsible for the satellite AGNs must 
occur when the primary component is still the central galaxy of 
its own DM halo. If the event is immediately followed by a halo 
merger, and the remnant galaxy becomes a satellite, it is classi- 
fied as an AGN satellite. Thanks to its implementation of strip- 
ping and the modeling of satellite-satellite mergers, a population 
of "true" AGN satellites is present in the KVB08 model. However, 
as shown in fig. [5] the contribution from this channel is quite small, 
i.e. f(AGN\P) are similar to those predicted by the other models. 

The KVB08, WDL08, SHC08 models do not allow "radio- 
mode" feedback in satellite galaxies by construction. This does not 
hold for MORGANA that provides a population of satellite radio 
galaxies with roughly the correct distributions. In addition, in this 
model f(LRA\M*) and f(HRA\M*) have similar shape and the 
fraction of satellites with AGN activity does not vary significantly 
as a function of halo mass. This is due to the assumed modelling of 
the cold gas reservoir around SMBHs: recently accreted satellites, 
despite being instantaneously stripped of their hot gas, are still able 
to sustain a short period of accretion, until their cold gas reservoir 
is completely depleted. 

In fig. [6] and [7] we show the bidimensional distributions 
f(A\MoM, M+) for central and satellite galaxies respectively. 
These figures provide a useful summary of the behaviors just dis- 
cussed. In these figures, darker grey levels correspond to increasing 
f(A\P), and are normalized to their highest value in each activity 
class. In the leftmost column, we show the observed distributions 
(as in P09, see their fig.6). In each model panel, we show the con- 
tours corresponding to the 25%, 50%, 75% number density lev- 
els in observed data. The strong correlation between stellar mass 
and halo mass for central galaxies is clearly visible in fig. [6] and it 
prevents us from decoupling the role of environment from that of 
stellar mass, as discussed in P09. Fig.[6]confirms that all models are 
able to reproduce the different distributions of SF and HRA centrals 
in the M* -Mum space: the former class is dominated by low-mass 
centrals, living in small DM haloes, while HRA galaxies are pref- 
erentially associated with more massive systems. AGN centrals in 
KVB08, WDL08 and MORGANA tend to populate the upper right 
region of the diagram, while AGN centrals in the SHC08 model 
are more likely associated with intermediate mass systems. Finally, 
LRA and HRA centrals populate different regions of this space in 
the models, in contrast with observational measurements. For all 
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contours in each model panel mark the 25%, (thin lines) 50% (solid lines) and 75% (thick lines) number density levels for the same activity class in P09 data. 



models, the overall agreement with the observed distributions is 
rather poor. 

For satellite galaxies, the activity seems to be more strongly 
related to stellar mass than halo mass, in qualitative agreement with 
observational data. In particular, fig. [7] highlights that active satel- 
lites are preferentially found in regions of the M* — Mdm diagram 
which are contiguous to those populated by active centrals belong- 
ing to the same activity class. This confirms that, in all models, 
active satellites are recently accreted objects. A different treatment 



for hot gas stripping in newly accreted satellites, as proposed in 
KVB08, leads to a better description of the distributions of SF satel- 
lites, which in this model populate a wider area in this diagram and 
show a stronger dependence on Kh than on A^dm- 



AGNinSAMs 



11 



© 



00 

o 



> 80% 

P09 



Satellites 

60%-80% ■ 40%-60% ■ 20%-40% < 20% 

WDL08 MORGANA KVB08 SHC08 



1— I 

o 

CV2 

o 

CV2 




CV? 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



Hf MM 1 1 1 1 +H 



III MM MM III 



LU 



1 1 1 1 ■ 1 1 1 1 1 1 



LLL 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




Willi MM M l 



III M ill M ill III 



in 



I ■ i ■ 1 1 1 1 1 1 1 



LLL 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



III MM IIII-H4 



III MM MM III 



LLL 



1 1 1 1 1 1 1 1 1 ■ I 



LLL 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 



f++ 1 1 1 1 MM III 



III M ill M ill III 



in 



1 1 1 1 1 1 ■ i ■ 1 1 



LLL 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



H+ 1 1 1 1 I III +K 



III M ill M ill III 



LLL 



ImiiIiimI 



ill 



12 13 14 12 13 14 12 13 14 12 13 14 12 13 14 

Log(M DM /h) (M ) 

Figure 7. Same as fig.[6]for satellite galaxies. Colours, shading contours and legend as in fig. [6] 



OS 

CO 



2 

o 
< 



< 



< 
OS 



5 SUMMARY AND DISCUSSION 

In this paper, we compare predictions from four different semi- 
analytic mo dels for galaxy forma tion with the observational data 
discussed in |Pasqualietal] ( l2009h . These data allow, for the first 
time, a careful comparison with theoretical models, also in terms 
of the galaxy hierarchy (i.e. their nature of centrals or satellites). 
This is important for a better understanding of the limitations and 
successes of current models, in which central and satellite galax- 
ies correspond to rather distinct evolutionary paths. The theoretical 
models considered in this study adopt different assumptions for the 
various modes of AGN activity: all of them distinguish between a 



bright phase (usually referred to as "bright-mode" or "QSO-mode") 
and a low accretion efficiency phase which is usually associated 
with the development of radio jets able to offset cooling flows 
(the so called "radio-mode"). In this paper, we focus on four ac- 
tivity classes, defined by the star formation rate, optical AGN ac- 
tivity, and radio emission associated with the radio-mode. Model 
accretion rates onto SMBHs are converted into 'observables' by 
means of empirical relations calibrated on local samples. We have 
attempted to reproduce the observational selections and checked 
the influence of different selection criteria: the absolute fractions 
of galaxies in each activity class as a function of both stellar and 
parent halo mass depend on the details of both the selection criteria 
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and the adopted conversions. The shapes of the predicted distribu- 
tions, however, are not significantly affected. 

We show that the four models used in this study predict in a 
few cases quite different distributions as a function of stellar and 
parent halo mass. This is interesting and demonstrates that the ob- 
servational measurements can discriminate between different phys- 
ical models. In particular, we have shown that in all models SFGs 
are typically low-mass galaxies residing in low-mass haloes, while 
powerful radio galaxies are expected to reside at the center of rel- 
atively massive haloes. Although this is in qualitative agreement 
with observations, we highlight a number of model "failures": 

• almost all massive centrals are predicted to host a detectable 
radio source, at variance with observations. 

• The f(LRA\P) and f(HRA\P) distributions differ from 
those observed for all models. In particular, it is possible to iden- 
tify a typical halo mass separating a 'low' from 'high' radio activ- 
ity (LRA and HRA classes). The precise value of this 'transition' 
halo mass depends on the adopted modelling for AGN feedback. 
We stress that in the WDL08 and KVB08 models this result is not 
surprising, since accretion rates onto BHs during radio-mode are 
explicitly related to the mass of the parent halo (see eq. [2] and [5] 
respectively). 

• The distributions of f(AGN\P) are generally too skewed to- 
wards high stellar and halo masses with respect to the observational 
measurements. This results in activity levels that are too high for 
massive central galaxies in clusters. The SHC08 model is an ex- 
ception as the predicted distribution in this case is skewed towards 
low-mass galaxies and low-mass haloes. In addition, this models 
predicts total fractions of galaxies with AGN activity that are signif- 
icantly lower than those predicted by the other three models used in 
this study. We ascribe these differences to the inclusion of a critical 
mass for BH accretion: since at low redshift most massive galaxies 
have central SMBHs with masses already larger than their critical 
mass, further accretion onto these objects is suppressed. 

• The distribution of star forming centrals in the Af* — Mom 
plane is well reproduced by the WDL08, KVB08 and SHC08 mod- 
els, while MORGANA predicts a wider distribution of SF galaxies, 
extending towards larger stellar and halo masses. This is due to an 
inefficient quenching of star formation in massive haloes. 

• The KVB08 model exhibits the best agreement with the ob- 
served distributions of star forming galaxies. In particular, the 
model predicts the correct dependence of the f(SFG\M±) on M* 
(although with a slightly higher normalization), both for centrals 
and satellites. As the other models used in this study, however, it 
does not reproduce the observed distribution of SF satellites as a 
function of parent halo mass. 

The disagreement between the predicted and observed frac- 
tions of radio sources represents a potentially serious problem for 
current models of galaxy formation and evolution. The inclusion of 
"radio-mode" AGN feedback in theoretical models of galaxy for- 
mation was originally motivated by the need for a powerful energy 
source, that would be able to counteract gas cooling at the centres 
of massive haloes. Most models assume a strong dependence of 
radio-mode feedback on the parent halo mass: this is essential in 
order to reproduce the observed rapid cutoff at the bright/massive 
end of luminosity /mass function, and the old stellar population ages 
of massive galaxies. As a consequence of this assumption, these 
models predict that essentially all massive galaxies should be asso- 
ciated with a bright radio source, while observational data suggest 
that faint and bright radio sources are found in similar environ- 
ments in equal numbers. If, as suggested by our analysis, we need 



to assume a weaker dependence on halo mass, then this physical 
process might not be capable of offsetting the cooling flows alone. 
This is hinted at by the fact that the model which predicts the low- 
est fraction of radio sources (MORGANA) is, at the same time, the 
model which provides the worst agreement with the distribution of 
star forming galaxies. 

Details of the adopted modelling for BH growth (e.g. the duty 
cycle) might also play a role. In particular, the models used here 
neglect completely one relevant aspect of AGN activity: the pres- 
ence of well defined duty cycles of radio sources (of the order of 
10 7 - 10 8 yrs for radio-loud galaxies). In SAMs, both the "radio- 
mode" feedback and the cooling rate are quantities integrated over 
a well defined time interval. The "quenching of the cooling flow" 
is therefore assumed to happen on this time-scale. T hese assump- 
tions are supported by results from lBest et al.l d2006h , who showed 
that the time-averaged energy output from recurrent radio-sources 
in elliptical galaxies balances almost perfectly the energy losses 
from the hot gas in the parent DM halo. If we assume that the radio 
sources are active only for a fraction of the adopted time-step, the 
number of sources detected as radio galaxies will be reduced. At 
2 ~ 0.1, however, the typical time-steps adopted in all SAMs we 
consider are shorter than or of the same order as the duty cycle of 
radio loud galaxies, so that the expected correction should be small. 
Even for models with integration time-scales larger than the typi- 
cal radio-loud duty cycle, it is not clear if a shorter "radio-mode" 
feedback activity would be as efficient in quenching cooling flows. 
In fact, the present approach assumes the heating term to perfectly 
offset the cooling when the "radio"-mode is "on", while the ther- 
mal evolution of halo gas proceeds as if there had never been a 
radio source present when the "radio"-mode is "off". This implies 
that introducing a duty cycle in these models would probably re- 
duce the efficiency of quenching cooling flowo m addition, we 
note that even fixing the total fractions of radio galaxies, we would 
still have the problem that their distributions do not reproduce ob- 
servational measurements. Therefore, more realistic scenarios have 
to be considered. These should take into account, for example, that 
radio jets may add entropy to the gas, and that after a complete 
radio cycle, the halo gas would be hotter and less dense than be- 
fore the onset of the "radio"-mode. Cooling will thus proceed at 
a slower rate. Our models also do not take into account any pos- 
sible "excess" in the radiated energy (with respect to the energy 
needed to quench cooling flows), which could unbind the gas from 
the halo, again sup pressing cooling on a longer timescale (see e.g., 
iBower et al .1120081 )- Therefore, an accurate modelling of AGN duty 
cycles requires a deep revision of AGN feedback schemes in SAMs 
that goes beyond the aims of this paper. 

Our results suggest that either 1) more sophisticated modelling 
of the triggering of radio activity and its impact on the thermo- 
dynamics of the surrounding gas is needed or 2) other processes, 



6 We note that MORGANA already implements some of these ideas (al- 
though in a very simplified way). This model explicitly follows the build up 
of a gas reservoir around SMBHs as a result of cold gas infall. Gas is then 
accreted onto the central object on a time-scale determined by the viscosity 
of the accretion disc, and originates collimated jets, that heat the halo gas. 
The time-interval between the onset of the cooling flow and its quenching 
through jet heating could be considered the equivalent of a 'duty cycle'. We 
have shown that this model performs slightly better in terms of the radio 
active fraction (not all massive galaxies host a bright radio source), but it 
is not able to quench cooling flows as efficiently as needed to reproduce 
the observed colours and star formatio n rates of massive galaxies (see e.g. 
iFontanot et alj200llKimm et alj2009h . 
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Figure 8. Fraction of central galaxies belonging to the activity classes HRA 
and LRA as a function of the mass of the parent DM halo, for different 
choices ofip = ikin/^RM- In a U panels lines and shading as in fig. [4] 



such as gravitational heating by infalling satellites or gas clumps 
dKhochfar & Ostrikenl2008l [pekel & Birnboimll2008l) . are impor- 
tant in preventing the gas in massive halos from cooling. 

In all the models, active satellites represent the tail of the cor- 
responding distributions computed for central galaxies. This is due 
to our simplified treatment of satellite evolution. In particular, three 
out of the four models used here assume that the hot halo associated 
with a galaxy is instantaneously stripped once the galaxy becomes 
a satellite. The cold gas available is rapidly turned into stars and 
ejected in the surrounding medium, never being reicorporated onto 
the same galaxy. The KVB08 model represents an improvement, as 
it allows satellite galaxies to keep their hot gas reservoir for some 
fixed time after accretion. We have shown, however, that this model 
does not reproduce the correct dependence of SF activity as a func- 
tion of halo mass. 

Our work confirms that statistical studies aimed at disentan- 
gling the role of stellar mass and parent halo mass (i.e. the "nature" 
versus "nurture") are hampered by the pivotal role played by central 
galaxies in the models. An improved description of the evolution 
of satellite galaxies is of fundamental importance to understand the 
relative importance of various physical processes responsible for 
galaxy activity in different environments. 

A number of additional caveats should be considered. In order 
to compare model predictions with observational data, we had to 
apply a number of corrections and empirical conversions. In partic- 
ular the adopted conversion from estimated accretion rates into ra- 
dio luminosities is critical. To test the robustness of our results, we 
repeated our analysis by assuming that only a fraction tp of the heat- 
ing energy is sufficient to originate jets and bubbles (parametrized 
as ip — I/kin/iRM). With respect to the results presented in fig.[3] 
this implies lower radio luminosities associated with the "radio- 
mode", but preserves the shape of the predicted distributions as 
a function of A/dm- We show in fig. [8] the f(HRA\MoM) and 



f(HRA\MuM) distributions corresponding to three different val- 
ues of ip. This figure shows that the number of sources in each class 
depends significantly on the value assumed for ip, but our main 
conclusions are unchanged: the distributions for LRA and HRA 
galaxies are different and there are too many HRA model galaxies 
with respect to observational measurements. It is however inter- 
esting that the peak of the LRA distribution tends to shift towards 
higher masses as ip decreases, thus slightly reducing the disagree- 
ment with the observed distributions. 

Our analysis highlights new problems for current models of 
galaxy formation and evolution, and in particular for currently 
adopted models of AGN feedback. Although this mechanism helps 
to solve a number of long-standing issues, it provides a poor match 
to the observed fractions and distributions of radio galaxies. From 
an observational point of view, a better and more complete sam- 
pling of the activity levels as a function of different physical prop- 
erties and environment, would be crucial to confirm and strengthen 
our conclusions. Important additional information will come from 
relating galaxy activity with other physical properties. For exam- 
ple, recent studies have suggested that the relation between ac- 
tivity and morphological type can provide interesting c onstraints 
on the accretion mechanisms dGeorgakakis et al.l 120090 , and on 
the co-evolution betwee n host galaxies and their central SMBHs 
dSchawinski et alj|201fj|) . 
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